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^ , A new general formalism for determining the electric multipole polarizabilities of 

quantum (atomic and nuclear) bound systems based on the use of the transition 
matrix in momentum space has been developed. As distinct from the conventional 
| approach with the application of the spectral expansion of the total Green's function, 

our approach does not require preliminary determination of the entire unperturbated 
spectrum; instead, it makes possible to calculate the polarizability of a few-body 
bound complex directly based on solving integral equations for the wave function of 
the ground bound state and the transition matrix at negative energy, both of them 
being real functions of momenta. A formula for the multipole polarizabilities of a 
two-body bound complex formed by a central interaction potential has been derived 
and studied. To test, the developed T-matrix formalism has been applied to the 
" calculation of the dipole, quadrupole and octupole polarizabilities of the hydrogen 

atom. 
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1. Introduction 

Study of the few-body quantum systems persists a perspective line of the advance- 
ment of the modern physics. In atomic physics, the Efimov effect for the three-body 
system, the existence of which was predicted in 1970 (Ref. 1), has been experimen- 
tally verified using ultracold atoms and tuning the atom-atom scattering length by 
the magnetic field near a Feshbach resonance 2 (see, for example, the reviews 3-5 ). In 
nuclear physics, extensive investigations of few-body nuclei advantageously under- 
taken over the past fifty years make it possible to gain new important information 
about the nuclear force. 

The investigation of the behaviour of the few-body nuclei in the external electro- 
magnetic field permits to obtain additional data on their properties, specifically, on 
the electric polarizabilities and magnetic susceptibilities as important fundamental 
quantities of complex systems. In spite of the fact that the study of the electric 
polarizabilities of the few-body nuclei has attracted considerable interest from both 
the experimental 6-11 and theoretical 11-25 points of view and definite progress in the 
field has been reached, a great deal needs to be done in this area. 

Up to now, there exists a discrepancy between the direct experimental result 
for the electric dipole polarizability of the nucleus 3 He (determined by measuring 
deviations from the Rutherford scattering law of the low-energy elastic 3 He 
scattering by the Coulomb field of the heavy nucleus 208 Pb (Ref. 8)) and the result 
deduced from the data for the total 3 He photoabsorption cross section using the sum 
rule a_ 2 (Ref. 9). 
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The conventional formula for the electric dipole polarizability of the iV-particle 
bound complex has the form 

«si = 2 2^ 7T— 5 > I 1 ) 

where Di = eii"i is the dipole moment operator, and ^ are the energy and the 
wave function of the z/-excited state, ri and p are the radius vectors of the particle 
1 relative to the center of mass of the bound complex and the center of mass of 
the complex relative to the charged particle creating the electric field, and the 
summation is taken over all possible excited discrete bound and continuum states. 

Although the formula (1) using the spectral expansion can be practically applied 
in the case of the two-body complex (N = 2), the direct determination of ckei 
for three and more body complexes (N > 3) by this relation with the performed 
beforehand calculation of all discrete and continuum states is not feasible. The 
known method of Dalgarno and Lewis 26 permits to bypass the difficulty related 
with taking account of the intermediate continuum states in all possible channels in 
Eq.(l) introducing an additional function | F), 

= 2(*„ | Dl - a, I , I ^ E 1 *" )( *" 1 *V° 1 *°> ■ 

where | F) being the first-order correction to the unperturbated state satisfies an 
inhomogeneous differential equation. Dalgarno-Lewis method has been effectively 
applied in Refs. 12 and 23 to calculate the electric polarizabilities of hydrogen and 
helium isotopes. Though the result for a_ei( 3 He) obtained therewith supports the 
value deduced using the sum rule <r_2 (Ref. 9), the subject of the discrepancy 
of the experimental data for 3 He under discussion may not be considered as 
conclusively established since indispensable consistent calculations of the electric 
dipole polarizabilities of the three-body nuclei on the basis of the rigorous 
mathematical Faddeev's 27 formalism have not be performed yet. 

In this paper we formulate the direct T-matrix approach (DTMA) to determi- 
nation of the electric polarizabilities of a bound system, that relies on the solution 
of the integral equations for both the bound-state wave function and corresponding 
components of the partial transition matrix of the system. In Section 2, following 
the Watson-Feshbach method, we express the polarization potential describing in- 
teraction between a charge particle and a bound complex consisting of N particles 
in terms of a "truncated" Green's operator of the system. In this way we obtain a 
general expression for the electric multipole polarizability of the system. In Section 
3 the T-matrix approach to determination of the electric polarizabilities is formu- 
lated and simplifications of the general formula assuming conservation of the space 
parity and the total orbital moment of the system, are considered. Section 4 is de- 
voted to the application of the elaborated formalism to the two-body bound systems 
with the central interaction. It is shown that the electric 2 A -pole polarizability of 
the two-particle 5-state bound complex contains information both on derivatives 
(of the order A and lower) of the wave function and on the partial component of 
the transition matrix that corresponds to the orbital state with I = A. Section 5 
contains the application of the T-matrix formalism to determination of the electric 
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multipole polarizabilities of the hydrogen atom, discussion and conclusions. 
2. Polarization potential 

The formula for the polarization potential describing the interaction between a 
charged particle and a bound complex consisting of N particles follows imediately 
when treating the N + 1 body problem, namely, the low-energy scattering of the 
complex by the Coulomb field of the charged particle with the kinetic energy of 
the relative motion of the particle and the complex e being well below the breakup 
threshold energy of the complex. Experimentally, to determine the polarizability of 
the complex in the direct way a heavy nucleus (the particle 0) is used as a source 
of the intense electric field. For simplicity sake assume that the iV-body complex 
contains only one charged particle 1, the other particles (2, 3, . . . N) being neutral. 
The Hamiltonian of whole system of N + 1 particle has the form 

H = h° + H + vg, (2) 

where is the operator of the kinetic energy of the relative motion of the particle 
and the center of mass of the complex, H is the total Hamiltonian of the iV-body 
complex, H = H° + V, H° is the kinetic energy operator of the relative motion of 
the particles inside the complex, V is the total interaction potential of the complex, 
V = v ij , v oi is the operator of the Coulomb interaction between the charge 
particles and 1 (with the charges e and ei), is the operator of the interaction 
between the particles of the complex i and j. 

Using the known Watson-Feshbach projection technique 2 ' 29 ' 30 with the projection 
operators P =| >< I an d Q = 1 — P, where \l/ is the wave function of the 
ground bound state of the complex of iV particles (1, 2, ... N) with the binding 
energy B (normalized to one), the operator of the effective interaction may be 
written as 

V eff (e) =< I n{E) | #o > , (3) 

where the averaging is taken over the variables of the relative motion of the particles 
of the complex. The operator 1Z{£) satisfies the integral equation 

n£) = vg + vgg9 t (£)n(£), (4) 

in which the potential energy operator is the "external" Coulomb interaction 
potential between the field source (the particle 0) and the charged particle 1 of 
the complex, vg, and the propagator is a "truncated" Green's operator containing 
the "internal" interaction potential of the complex V, 

<3?nt{£) = QQint(S), Q mt (E) = {E-W- V)~\ (5) 

where W is the free Hamiltonian of the complete system (of N + 1 particles), 
%° = h° + H°, S is the total energy of the system, £ = e — Bq. 

Introducing the transition operator T{£) which satisfies the integral Lippmann- 
Schwinger equation 

T{£) = V + Vg°(£)T(£), (6) 
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the operator Gint(£) can be written in the form 

G int (£) = g°{£) + g°(£)T(£)g°(£), (7) 

where the free propagator G° is given by G°(£) = (£ — W)^ 1 . 

The polarization potential that corresponds to the second order of the perturba- 
tion expansion (in powers of the Coulomb interaction vg) of the operator 1Z(£) in 
the expression for the effective potential (3) is given by 

V pol (e) =< * I vgG? nt (£)vg | * >, (8) 

where 

G? nt {£) = G°(£) + G°{£)T{£)G°{£) - PG int {£). (9) 

According to the uncertainty principle, in the case of asymptotically large (in 
comparison with the complex size) distances between the particle and the complex, 
Po, the momentum variables of relative motion among particles inside the complex 
considerably exceed the momentum variable of relative motion of the particle 
and the center of mass of the complex. In such a case, in the expression for the 
"truncated" Green's function £?£ t (£) in Eq.(8) we can neglect with the variable 
quantity that corresponds to the operator e — h° in comparison to the variable 
quantity that corresponds to the operator H° + B . Then the polarization potential 
takes the known local form 

2 oo 

V P ol( P 0) = ~% £ 4*Rf , (10) 
Z A=l PO 

where each of coefficients in the sum (10), C(e\, that characterizes the strength of 
the individual constituent of the polarization potential with the asymptotic Po 2A ~ 2 
is the electric polarizability of the iV-particle complex of the multipolarity 2 A 

a EX = -2 < | M X G Q (-B )M X | tt > ■ (H) 
Here M x is the multipole moment operator of the charged constituent particle 1, 

M x = e 1 r x 1 P x (r 1 -p ), (12) 

i - ! is the radius vector of the particle 1 relative to the center of mass of the bound 
complex; the radius vector p describes the position of the center of mass of the 
bound complex relative to the charged particle creating the electric field, P\(x) is 
Legendre polinomial (the unit vectors are marked with a hat, a = a/ a); the quantity 
G®(—B ) is the "truncated" Green's operator of the complex G Q (E) = QG(E) at 
the bound state energy E = —B , G(E) = (E — H° — V)^ 1 is the total Green's 
function of the complex. 

Note that the formula (1) readily follows from the general formula (11) after 
applying the spectral expansion of the total Green's operator G(—Bq) in the 
complete set of the eigenfunctions of the total iV-particle Hamiltonian H. 
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3. T-matrix method of determination of the electric multipole 
polarizabilities of a bound complex 

In this Section we formulate a method for the direct calculation of the polarizability 
of a few-particle bound system starting immediately from the definition (11) without 
recourse to expansion in terms of a set of excited discrete and continuum wave 
functions. 

The general formula (11) for the electric dipole polarizability of a bound 
complex of iV particles (1, 2, ... N) contains the ground-state wave function and 
"truncated" propagator 

G Q {E) = G(E) - PG{E) (13) 

at the energy of the bound state E = —B . The wave function \l/ satisfies the 
Schrodinger equation 

(H + V + B )^ = 
or the equivalent homogeneous integral equation 

= G°(-B )V* . (14) 

We shall restrict our consideration to the case that the complex is in the ground 
bound state. 

The "truncated" propagator G Q (E), as distinct from the total one, G(E) does 
not contain the pole singularity at E — —B . In an explicit form, we find the 
expression for G Q (E) writing the total Green operator G(E) in terms of the free 
propagator G°(E) = (E — H°)~ l and the transition operator T(E), 

G(E) = G°(E) + G°(E)T(E)G°(E), (15) 

where the operator T(E) is determined by the Lippman-Schwinger integral equation 

T(E) = V + VG°(E)T(E) . (16) 

In view of Eq.(15) and the equality 

pg{e) = eTb/' (17) 

the expression (13) takes the form 

G Q (E) = G°{E) + G°{E)T(E)G°{E) - - 1 P . (18) 

E + B 

Separating out from the total transition operator T(E) the singular pole part 
that corresponds to the ground bound state of the system at the energy E = —B , 
we write 

T(E) = 1 ' + T(E) , (19) 

where the vertex function | T ) is expressed through the wave function of the ground 
bound state of the system, 

| r > = [G^-Bo)]- 1 I * > = V | *o> , (20) 
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satisfying the homogeneous integral equation 

| r > = vg°(-b ) I r > , (21) 

and T(E) denotes the smooth (non-singular) part of the transition operator. 

Taking into account the expression (19) we write the formula (18) for the operator 
G Q {E)m the form 

Gfl{E) = G°(E) + G°(E) [G ° { ~ B ° )rl ' ^°| ( ^ ' [GO{ ~ BoTl G°(E) 

E + B + G °( E )T(E)G°(E) . (22) 

Cancelling the pole terms in Eq.(22) with the use of the identity 

G°(E) \G\-B )}- 1 = [G^-Bo)]- 1 G°(E) = 1 - (E + B )G°(E) (23) 

we may write the expression (22) in the form 

G Q (E) = G°(E)-G°(E)\y )(y \-\y )(y \G°(E) 

+ (E + B )G°(E) | ^ )(^o I G°(E) + G°{E)f(E)G°{E) . 



(24) 



At the point of the negative energy of the bound state of the A^-particle complex, 
E = —B , the operator G Q (—B ), that is contained in the general formula for the 
electric dipole polarizability (11), becomes 

GQ(-Bo) = G°(-B ) G (-B )\* )(* \-\* )(*o\G°(-B ) 

+ G°(-B )f(-B )G°(-B ) . (25) 

With the use of the expression (25) the formula for the polarizability (11) takes 
the form 

a EX = - 2(^ | M X G°(-B )M X | *„) 

+ 2(^ | M X G°(-B ) | ^ )(^o I M x | * ) 
+ 2(^ | M x | ^ )(^o I G°(-B )M X | * > 
- 2(^ | M X G°(-B )f(-B )G°(-B )M X \ * ) ■ 



(26) 



A noticeable simplification of the general formula for the electric multipole 
polarizability (26) takes place, if the interaction potential V is invariant relative 
to the space reflaction (for example, for the systems with the Coulomb or nuclear 
interactions). In this case the wave function of the bound complex is characterized 
by a definite parity. The conservation of the parity leads to nullification of the 
matrix elements 

(*o I M x | * > , (# | G°(-B )M X | * ) and (* I M X G°(-B ) \ V Q ) , (27) 

which are present in the second and third summands of the formula (26), at odd 
values of A. 



- 7- 



The formula (26) is also simplified in the case of the invariance of the interaction 
relative to rotations that leads to conservation of the angular momentum. In such a 
situation the matrix elements (27) are proportional to the Clebsh-Gordan coefficient 
(LXM/i | LM'), satisfying the triangle condition A(LXL). Specifically, they vanish 
after integrating in angular variables for all A ^ 0, if the total orbital moment of 
the complex is equal to zero, L = 0. 

In the case when the ground bound state of the complex is characterized with the 
total orbital momentum L = 0, that is real for the symplest atomic and few-hadron 
systems, the formula for the electric multipole polarizability (26) is simplified to the 
expression 

a EX = - 2(^ | M X G°(-B )M X | tf > 

(28) 

- 2(^ | M X G°(-B )T(-B )G°(-B )M X | tf > • 

As is evident from (26) or (28), to determine the electric dipole polarizability 
of the bound complex, it is necessary to know not only its wave function, but 
the smooth part of the transition matrix at negative energy of the bound state 
(E = —B ) as well. The corresponding transition matrix can be determined 
by solving the Lippmann-Schwinger integral equation for the two-particle system, 
the Faddeev integral equations 27 for the three-particle system or the Faddeev- 
Yakubovsky equations 28 for more complex systems. 



4. Electric multipole polarizabilities of the two-particle bound complex 

In the case of two-particle bound complexes, the derived formula for the electric 
multipole polarizability (28) is simplified. Considered here is a stable bound complex 
consisting of a charged particle 1 and a neutral particle 2. The interaction potential 
between the particles is taken to be central, the complex is in S'-wave ground 
bound state. We denote quantities describing two-particle system by small letters 
as distinct to notations by capital letters of quantities for iV-particle complexes 
used in the foregoing Sections: b = ^ 2 Kq/2jUi 2 and -^o are the binding energy 
and the wave function of the two-particle bound complex, g°(—bo) = {—bo — h )^ 1 
is the free Green's operator, t(—b ) is the two-particle transition operator, /i 12 = 
mim2/ '(mi +m 2 ) is the reduced mass of the particle 1 and 2. Further consideration 
we perform in the momentum space. 

Starting from (28) and taking into consideration that the smooth part of the 
two-particle transition matrix t(— bo) has form of the sum of the smooth part of its 
S'-wave partial component (with / = 0), to(-bo), and the sum of all higher partial 
orbital components (with / > 0), t^(—b ), which are non-singular at the energy of 
the bound state, 

t(-bo)=t (-bo)+t w (-bo), (29) 

we write the formula for the electric multipole polarizability of the two-particle 
system in the form 

a EX = - 2(^ | M x g°{-b )M x | ^o) 

( 3 °) 

- 2(^o I M x g°(-b ) t (-bo)+^ h) (-bo) g°(-b )M x | Vo> • 
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The wave function of the complex in the S'-wave bound state satisfies the 
homogeneous integral equation 

I ^o) = 9°(-b )v | ^o> ■ (31) 

and each of the partial components of two transition matrix ti(e) satisfies the 
inhomogeneous integral Lippmann- Schwinger equations 

ti(e) =vi + v l g°(e)t l (e) . (32) 

where vi is the partial component of the interaction potential, e is the energy of the 
relative motion of the particles. 

The solution of the equation for the partial transition matrix (32) that 
corresponds to the orbital moment of the ground state ((/ = 0)) has the form of 
the sum of the pole and smooth operators 

*o(e)= l7 °^; 7ol +to(£), (33) 
e + bo 

where the vertex function | 70) = i>o | ipo) = ( — ^0 _ h°) \ ipo) satisfies the 
homogeneous integral equations that follows from the equation (31), 

I To) = v g o (-b ) I 70) . (34) 

Using the operator wo( e ) = ( e + bo)t (e), which satisfies the inhomogeneous 
equation 

w (e) = (e + b )v + v g°(e)w (e) • (35) 
we write the smooth part of the transition operator, to(e) in the form 

7 f x ^o(e)- I 7o)(7o I 

^(e) = —t • (36) 

e + Oq 

At the point e = — 6 , the inhomogeneous equation (35)becomes homogeneous one, 

wo(-6o) = v g°(-b )cj (-b ) . (37) 

Since the kernels of the equations for the function | 70) (34) and the operator u (— 6 ) 
(35) coincide, the solution of the operator equation (37) may be written as 

wo(-6o) =| 7o)(7o I • (38) 

According to (36), the operator t (— b ), which is contained in the formula for 
polarizability (30), is related to the operator cuo(e) by the expression 



7 r t \ duo(e) 
to{-bo) = 



de 



= u'oi-bo) . (39) 

-bo 



The operator u' (e) can be deduced by performing the differentiation of the equation 
(35) with respect to e and the inverted transition from v to t and Uo, 

^( e ) = " o(e)+ " o(e) f (£)a;o(e) , g^e) = -(e-hT 2 . (40) 
e + b 
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Evaluating on the right-hand side of (40) the inderterminacy of the type at the 
point e = — b Q , which appears in view of the relation (38) and the normalization 
condition (70 | g°\— bo)] | 70) = — (ipo \ ipo) = —1, we obtain the expression for the 
operator t ( — bo), 

^(-fco) = 4(-6„) =| 7o> (~) <7b I • (41) 

where 

Ri = -^o(7o I 9°"(-bo) I 7o) , 9°"(e)] = 2( e - h°y 3 ■ (42) 

In the case of the central interaction, the part of the transition matrix t^ h \—bo) in 
(30) may be written in the form of expansion in the set of the spherical functions of 
the angular variable momenta (with / > 1), 

00 1 

< k|^)(-& )|k' >=4vr£ ]T t l {k,k'--b )Y lm ^)Y^) . (43) 

1=1 m=—l 

Substituting the expressions (41) and (43) for t (-b ) and t h (-b ) into the 
formula for the electric polarizability (30) and taking account of the action of the 
multipole moment operator M x (11) on the function ifj Q in the momentum space 
(with r = iWk), 

(k I M x I ^ ) = i x e, (-H±\\ x (k)P x (k-p ) , (44) 
V mi2/ 

where 

<p x (k) = (-l) x k 

we perform integration with respect to the angular variables. 

Notice that the contribution of the smooth part of the S'-wave components of 
the transition matrix, t (— b ), which is of the separable form (41), proves to be 
equal to zero as a result of the conservation of the space parity and the orbital 
moment of the relative motion (similar to zero contributions in the general case 
from other factorable terms of the "truncated" Green's operator (25), see comments 
below the formula (26)). Non vanishing contribution in (30) makes only the partial 
component with / = A from Eq.(43), t\(k, k';— bo) that satisfies the Lippmann- 
Schwinger integral equation (32), 

/•oo dk"k" 2 1 
t x (k, hf- -b ) = v x (k, k') - / ^^vx(k, k")— 2 —t x (k", k'; -b ) , (46) 

where v x (k, k') is partial component of the interaction potential. 

The final result for the electric multipole polarizability of the two-particle 
complex with a central interaction between particles following from Eq. (30) is 
written in the form 



[Ilk) *><*> 



77712 = 777i + 777 2 , 



(45) 
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(2A + 1)tt 2 h 2 \m 12 J [Jo k 2 + n 2 

1 m i m 2 j hh 2 [°° AW2 <Pl(k)tx(K k'; -b )ip x (k' 
n 2 h 2 m 12 Jo Jo (k 2 + K 2 )(k' 2 + «g) 

According to (45), the functions <p x (k) in (47) are expressed in terms of the 
derivatives of the wave function of the ground bound state ipo(k) with respect to 
the relative momentum variable k. In the cases of the dipole (A = 1), quadrupole 
(A = 2) and octupole (A = 3) polarizabilities they are of the form 



<P\{k) 



-il>' (k) for A = 1 , 

+<(k) ~ Wo{k) for A = 2 , (48) 

+ iKik) - ±%{k) for A = 3 . 



The formula (47) derived in the case of the central interaction between the 
constituents of the two-body complex demonstrates that the electric 2 A -pole 
polarizability of the complex contains information not only on the derivatives (of 
the order A and lower) of its wave function, but on the partial component of the 
transition matrix, t x (k, k'; —b Q ), as well. 

It is advantageous to write the formula (47) in the form more convenient for the 
practical use. In order to do this, instead of the two-particle t-matrix (that is the 
function of two momentum variables k and k'), we introduce in (47) the function of 
one variable 

roo dk'k' 2 1 
Uk) = ~ / -i^t x (k'k'; -b )— 2 -<p x {h>) (49) 

that satisfies the inhomogeneous integral equation with the kernel of the Lippmann- 
Schwinger equation (46) 

roo (jk'k 12 1 

Mk) = fx(k) - / ~^v x (k, k')— 2 — ——(f>x(k') (50) 

in which the free term is determined by the formula 

roo dk'k' 2 1 
h(k) = - ^v^kX)^^—-^). (51) 
Jo 2n 2 i^ 2 +bo 



In such a case, the formula (47) takes the form 



.hss! f^) 2A+1 r dke mkm^m . (52) 

\m 12 J Jo k 2 + k£ 



(2A+1)tt 2 h 2 \m 12 J Jo k 2 + Kq 



5. Application of the DTMA, discussion and conclusions 

Advantages of the t-matrix approach are manifested in calculation of polarizabilities 
of a quantum system, even if it is two-particle (the deuteron nucleus or the hydrogen 
atom). In the special case that the interaction between particles is purely «S-wave 
(for example, the separable interaction potential 31 ), in the formula (47) only the 
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first term persists. The corresponding reduced formula for the dipole polarizability 
of the two-particle bound complex has been earlier derived in the framework of 
the three-body formalism of the effective interaction of a charged particle and a 
complex 13 ' 14,32 . Applications of the t- matrix approach to calculation of the electric 
dipole polarizabilities of the deuteron, the triton and the lambda hypertriton as 
two cluster systems with purely »S-wave interaction have been carried out in the 
preceding our papers 24 ' 25 . 

For the deuteron with noncentral interaction between the proton and the neutron 
(in the state with the total angular momentum 1) the polarization in the electric 
field is anisotropic. Using the separable tensor potential 33 the longitudinal and 
transverse (relatively to the direction of the electric field) components of the dipole 
polarizability of the deuteron have been calculated in the work 15 . The results of the 
further calculations of the components of the deuteron electric dipole polarizability 
obtained in the framework of the chiral effective field theory 18 are in agreement with 
the results of Ref. 15. 

In the case of the hydrogen atom H (assuming that the proton mass is infinitely 
great when compared to the electron mass), the exact values of the electric multipole 
polarizabilities a E \(R) are known 26 . This is why it is possible to test directly the 
validity of the general formula for the polarizability of the two-particle system 
(47) - - both of the term with free propagation in intermediate state and of 
the terms with the multiple scattering in the higher orbital intermediate states 
(with / = 1,2,3,... for the dipole, quadrupole, octupole and higher multipole 
polarizabilities, respectively). 

Inserting the analytical expression for the partial components of the Coulomb 
transition matrix, obtained from the representation for the three-dimensional t- 
matrix < k|t c (e) |k' > derived in Ref. 34 with the use of the 0(4) rotation symmetry 
in four-dimensional Fock space 35 , into the formula (47) and separating out the Born 
term from the Coulomb t-matrix, 

t\(k, fc'; -6q) = v°(k, k') + Atf (*, k'\ -b ) , (53) 

we write o:e\(H) as a sum of three terms, 

a EX (R)=a° x + af + a^ s , (54) 

where a° x is the first term in Eq.(47) that describes free virtual propagation, is 
a part of the second term that contains only the Born term of the partial Coulomb 
t-matrix, v%, describing the single scattering in the intermediate state, and a^ s 
takes into account the multiple scattering contributions (of the order 2 and more). 

We show in the Table 1 the values of of the components a° x , af , ajf s and 
their sum (54) (in a. u.) that determines the electric dipole (A = 1), quadrupole 
(A = 2), and octupole (A = 3) polarizabilities of the hydrogen atom, a SA (H), 
obtained applying the direct t-matrix approach, together with the exact values of 
the polarizabilities derived by Dalgarno and Lewis 26 . Here, values of the quantities 
a\ and af are derived from (47) analytically, and the values of a{^ s are obtained 
numerically calculating the integral with At x (k, k'; —b ) that describes the multiple 
scattering in intermediate states. 
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Table 1. The components af and a^ s determining the electric dipole (A = 1), 
quadrupole (A = 2) and octupole (A = 3) polarizabilities of the hydrogen atom, aE\(K), calculated 
with the useof DTMA (Eqs.(47), (53) and (54)) together with the exact values ofthe polarizabilities 
(taken from Ref. 26) (in a.u.) 



A 


ofi x 






c*£;a(H) 


Exact 
value 
(Ref. 26) 


1 


7 

3 


10 

9 


1.055555... 


4.499999 . . . 


9 
2 


2 


81 

8 


105 
32 


1.593750 


15.000000 


15 


3 


99 


243 
10 


7.950000 


131.250000 


525 
4 



Data given in Table 1 indicate that the dipole polarizability of the hydrogen atom 
consists of two nearly equal parts — the first term with the free propagator, oi^ — X, 
and the sum of two others with the P-wave component of the Coulomb transition 
matrix, af + a^ rs = ^P. For polarizabilities of higher polarity the contribution of the 
term with the free propagator a° x still further increases reaching 67, 5% in the case 
A = 2 and 75,4% in the case A = 3. The contribution of the Born term a x to the 
polarizability «^a(H) diminishes slowly with A being equal to 24.7% for A = 1, 21.9% 
for A = 2 and 18.5% for A = 3. The contribution of the term describing multiple 
virtual scattering, ajf s , lowers with increasing A accounting for 23.5% when A = 1, 
10.6% when A = 2 and 6.1% when A = 3. It is worthy of note that contribution 
of the sum of the analytically tractable terms, a° + af , to the polarizabilities of 
the hydrogen atom aE\(H) is prevailing, it accounts for 76.5% (if A = 1), 89.4% (if 
A = 2) and 93.9% (if A = 3) of the total amount. 

The elaborated approach to determination of the electric multipole polarizabili- 
ties that relies on the two-particle t-matrix is found to be rather effective compared 
to the conventional method that uses the spectral expansion of the Green's function 
of the system. As evidenced by the results of the spectral-expansion calculations 
by Castillejo at al 36 , 65.8% of the magnitude of the electric dipole polarizability of 
the hydrogen atom comes from the P-wave excited bound states. Inclusion of the 
all excited bound states account for 81.4% . The rest 18.6% is provided taking into 
consideration the continuum states. 

To summarize, as one of advantages of the approach developed here is that there 
is no need to calculate the individual contributions from all discrete and continuum) 
excited states. Instead, the determination of one or several (depending on the in- 
teraction mode) partial components of the transition matrix is now required. It is 
also essential that the transition matrix in the main formula (28) depends on the 
negative energy of the bound state being therefore a real function of momenta. The 
developed method can be immediately extended to more complicated interactions 
between constituents, specifically, to the case of tensor interactions that gives rise 
anisotropic polarization properties of the system. Finally, it is quite important that 
the proposed method is suitable for more complicated three- and iV-body systems 
described by the Faddeev and Faddeev-Yakubovsky integral equations the kernels of 
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which are expressed through components of the transition matrices for subsystems. 
In the first place, we plan to apply the T-matrix approach to study the deformation 
properties of the 3 H and 3 He nuclei in the electric field. 
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